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Abstract 

For First Passage Percolation in Z d with large d, we construct a path connecting the 
origin to {x\ = 1}, whose passage time has optimal order log d/d. Besides, an improved 
lower bound for the "diagonal" speed of the cluster combined with a result by Dhar 
(1988) shows that the limiting shape in FPP with exponential passage times (and thus 
that of Eden model) is not the euclidian ball in dimension larger than 35. 
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1 Introduction 

Let {t(x, y), (x, y) edges of Z d } be a family of i.i.d. positive random variables. The quantity 
t(x, y) is called the passage time through edge (x, y). For a path V : xq —¥ x\ x n of 

neighbouring vertices, we denote by r(V) the passage time along V: t(V) = Y^i=i r ( x «-i; 
The family {r(x, y)} defines a random distance over Z d as follows: 

T>(x, y) = inf {r(V) ; V goes from x to y. } . 

We also set B t := {x; T>(0, x) < t}. 

This model is called first passage percolation. We refer to Kesten's St-Flour Lecture 
Notes [5] for a nice introduction to the subject. In this work we focus on the case where 
the common distribution of the passage times is the exponential distribution with parameter 
one. This case has received a particular attention for at least two reasons: 

• The process t i— > B t is then a Markov process. This is a consequence of the memory- 
lessness property of the exponential random variable. 

• Consider the random process B t only at random times at which a new vertex is added. 
Namely, set t = and t& +1 = inf {t > tk ', B t ^ B t _} and look at the sequence 
B to , B tl , . . . . This discrete process is known as the Eden growth process [I], which 
was introduced as a (very) simplified model for cells spread: to B tk one adds a new 
vertex x adjacent to B tk , with a probability which is proportional to the number of 
edges between x and B tk . 



We denote the canonical basis of M d by (ei, . . . , e^). For a subset A C M. d , the random variable 
T(A) is the first time the cluster B t hits A, that is, 

T(A) = min {t > ; B t n A ^ 0} . 

We will consider T(A) for the particular set H n = {x% = n}. Subadditivity theory was 
historically introduced by Hammersley and Welsh to show that the sequence T(n, 0, . . . , 0)/n 
converges almost surely to a constant \x = which is called the time constant. It can be 
seen as a consequence of the work of Cox and Durrett pQ that // is also the limit 

u = hm . 

a.s. n 

Kesten ([5], Th.8.2) was the first to prove that, for a large class of distributions for the r's, 
the constant /x is of order log d/d. His result was later improved by Dhar for exponential 
passage times. 

Theorem (Dhar [3]). For exponential passage times, 

lim u(g0 - - = -. 

Some numerical computations are included in [3], showing that the convergence is quite 
fast. Our aim in the present paper is to provide for a constructive proof of the fact that 
fj, = 0(logd/d), by exhibiting a path going from the origin to H% with small passage time. 
This path, except its last edge, is included in H . 



2 A path V connecting to {x\ = 1} 
A tree-like construction of paths 

Fix an integer i > 2, and pi,...,p£ a collection of i integers larger than 1. We want to 
associate to these integers a path V(£,pi, ■■■,Pe), having £+1 edges and connecting the origin 
to Hi, along which passage times are as small as possible. 

The first £ edges of this path lie in hyperplane H , and the last one connects H to Hi. 
The path V(£,pi, ■■■,pe) is defined as follows: 

Step 1 Among the 2(d — 1) edges (0, ± eA (with 2 < i < d), we consider the ones with the pi 
smallest passage times. This gives pi first edges ending at some vertices that we denote 
by x[l],...,x\pi]. 

Step 2 From each point £i[l], . . . , £i[pi], consider the 2(d — 1) — 2pi edges that are not collinear 
with one of the edges that have been already used in Step 1. Among these pi x (2(d — 
1) — 2pi) distinct edges we choose the P2 smallest ones. They end at some distinct 
vertices that we denote by ^[l], • • • , a^bs]- 
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Step £ From each point . . . , x^_i[p^_i], consider the 2{d — 1) — 2p 1 — 2p 2 — ■ ■ ■ — 2p e _i 

edges that are not collinear with one of the edges that have been already used in the 
previous steps. Among these pe-i x (2(d — 1) — 2pi — 2p 2 — ■ ■ ■ — 2pe_i) distinct edges 
we choose the pi smallest ones. They end at some distinct vertices that we denote by 
xt[l], ■ ■ ■ ,x t \pt]. 

Step £ + 1. Among the pi edges (a^[l], xe[l] + eo), • • • , (^[p^], xi\pi] + eo) we choose the one with 
the shortest passage time. We denote this path by xg — > X£ + \. 

Backtracking from xg + i to defines our path: this is the only path 

V(£,P!, ...,pt) : -)> x 1 ->■ x 2 -> x 3 ->■ ■ ■ • ->■ ->■ -)■ z m 

for which, for all 1 < % < £ — 1, the edge (xi,Xi+i) is of the type (xi[r], Xj + i[s]) for some 
integers r, s. 

The passage time of V 

Now our main result states that, for large d, one can find among the paths V(£,pi, ...,pe) a 
path whose passage time is of optimal order: 

Theorem 1. 

d e 
limsup- — - inf E[r(V{£,pi, -,Pi))] < t- 
d^oo logdi, P i,-,pt 4 

Proof. Let us first introduce a notation. Fix two positive integers n > k > and take 
e 1; e 2 , . . . , e n a family of i.i.d. exponential random variables with parameter one. Pick uni- 
formly one of the k smallest, we denote by f(n, k) its expectation. The mean passage time 
of V(£,pi, ...,Pi) can be written 

E[r(V(£,pi, ...,Pt))] = fW - l), Pl ) + f (px(2(d - 1) - 2 Pl ),p 2 ) + ... 

+ /L_ 2 (2(d-l)-2^ Pl ),p,J +f( Pi ,l). (1) 
From the well-known representation of the order statistics of ei, . . . , e n 

/ v (law) / 1 1 

(e (1) ,e (2) ,...,e (n) j = I -ei, -■ 
one readily deduces that 



1 1 1 

-ei H -e 2 , • • • , -ei H -e 2 H h e n 

n n — 1 n n — 1 



' ' 1 fc + 1 



^— ' ^— ' n — 7 2(n — k) 

i=l j=0 J v ; 



which implies from ([I]) that 



E[r(P(£, Pl , ...,p<))] < i(_^l±i— + tt^tUt n + 

2 v 2(rf-l)-pi 2(d-l)p 1 -(2pf+p 2 ) 



+ fLtl ) + I 

2(d - l)pe-i - (2p e _ 1 J2i<t-i Pi + Pi) Pi 
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We asymptotically minimize the above right-hand side by introducing a positive integer A 
and taking £ = [log d\ — A, and pt = [e l \ . We obtain 

l imsup ^L inf E[r(V(£, Pl ,...,pe))} < j [z V^a) 

d ^ +OQ \ogdt,Pi,..., Pe 4 yl-— e A J 

which gives the desired bound when A grows to infinity. □ 

If we use the same procedure to build a path from Hi to H 2 , H 2 to H 3 ,..., H n _i to H n , 
the family of passage times of these n paths is i.i.d. by construction. It then follows from 
the law of large numbers that 

limsup a- < e/4 = 0.679... 

d^oo log d 



Comments on the result 

1. We obtain a short and constructive proof of the bound \i < c st log d/d. We are however 
not able to achieve Dhar's optimal bound with the constant 1/2. The latter was 
obtained with a recursive argument applied to B t , but cannot provide for an effective 
path going to H n . 

2. Kesten's original proof of the existence of a path whose time constant is less than some 
constant (11 in his proof) times was also non constructive. However, it indicates 
that a path of length log d achieves this optimal order. This coincides with our choice 
of £ in the proof of Theorem 1. If one restricts the scope to paths of length 3 (£ = 2), 
one already gets an interesting bound i.e. 

limsuprfl inf E[r(7 3 (2,pi,p 2 ))] < C 

which proves that, in dimension large enough, the horizontal speed is bigger than the 
diagonal speed, which has been proved to be of order y/d (see next section). More 

£ 

generally, optimizing K[r(V(£,pi, ...,pe))] for a fixed £ leads to a bound n(d) < C /d 7 ^ 1 . 

3. Our result could be extended to a large class of distributions over the passage times, 
provided one has a good upper bound for f(n,k). This can be done if r has a first 
moment and a nice density near zero, not null at zero (such assumptions on r were 
considered by Kesten). 



3 Discussion on the diagonal speed and the limiting 
shape 

Richardson [7] proved that B t grows linearly and has a limit shape: there exists a nonrandom 
set B C R d+1 such that, for all e > 0, 

f((1-s)B C^c(1 + s)B ) 
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(Richardson didn't exactly deal with first passage percolation but with a class of discrete 
growth processes, including Eden's growth process.) The convergence also holds almost 
surely, the most general result is due to Cox and Durrett [1] . 

The shape Bo appears to be a ball of a certain norm, which is not explicit. Eden and 
Richardson observed on simulations that B looks circular in dimension two (though, they 
only performed the simulations up to a few hundreds vertices in B t ). Kesten has shown that, 
surprisingly enough, this is not the case for FPP with exponential passage times, at least 
when d > 650000. We conclude this paper by short arguments showing that B is not the 
euclidian ball when d > 35. Let us denote by /x* the "diagonal" time constant 

* y T(Jn) 

H* = hm 



a.s. n ' 



where J n = |xi + x 2 + ■ ■ ■ + Xd = ny/dj. Kesten ([5 J Th.8.3) observed that /i* > l/2eVd = 

0.184. ../y/d which, compared to /i — C(log d/d), gives that /x < /i* if d is large enough, 
yielding that B is not the Euclidean ball. Carrying his argument a little further, we obtain 
a slightly improved bound for fi*: 

Theorem 2. For all d > 2, 



> y/^T w 0,3313... 
~ 2Vd Vd 
where a* is the non null solution o/cotha = a. 

Proof. The proof is elementary, it mainly consists in bounding the probability that a fixed 
path going from the origin to J n has small passage times. 

We denote by the set of self-avoiding paths of length k (in the sense that they do 
not run twice through the same edge) starting from zero and hitting J n for the first time at 
time k. Because of self-avoidingness, passage times are independent along such a path. Fix 
a real number x > 0, since a path from the origin to J n has at least ny/d edges, 

P(T(j7n) < nx) < P( there exists k and a path V in s.t. r(V) < nx) 

= D k ] xP(r(M) <nx), (3) 

k>n\/d 

where T(k, 1) is a Gamma(/c, 1) random variable. The following estimate is straightforward: 

P(r(fc,l) <a)< [ae/kf . 
Lemma 3. For k ~ an\fd with some constant a > 1, 



card(C<"») (2d)* J ^ ( {a + 1)(M)/ ^ a _ 



1 We do not use the convention of Kesten [5] for the definition of /it*, yielding to a different factor of sfd 
between his statement and ours. 
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Proof of Lemma 0. We evaluate the number of (non necessarily self-avoiding) paths of length 
k which start from the origin and hit J n for the first time at time k. Such a path 

S : ->■ S x -> > S k 

is seen as one sample of the standard symmetric random walk in Z d . Its projection X = 
(0, Xi, . . . , Xk) on the axis x\ = X2 = ■ ■ ■ = Xd is a symmetric one- dimensional random walk 
with increments ±l/vd- Applying Th.9.1 in [6] and the Stirling formula gives 

(2dY k D ( ^ = P(S hits J n for the first time at k) 



J~( a - V 



□ 



Going back to (J3j) gives 



(nxe/k) f 



k>n\/d 

2\fdxe 



£ 2 E 

k>nVd 



27rnVd 1 (a + l)( a+1 )/ 2 «(a - 



where a = k/ny/d. This sum decays exponentially provided that 

x < sup {(a + l)( Q+1 )/ 2a (« - l)^- 1 )/ 20 } . 

2ey d a>\ 

This supremum is attained for the unique non null solution a* of cotha = a. It is equal to 
e-y/a 2 — 1. Theorem [2] then follows from the Borel-Cantelli Lemma. □ 

Corollary 4. In dimension d > 35, i/ie limiting shape is not the Euclidean ball. 

Proof. Combining Theorem [2] with Dhar's numerical computations we obtain for d = 35 that 
Ai(35) < 0.93 log(2d)/2d < 0.3313/v^ < /x*(35). □ 
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